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Improving quantitative dosimetry in ^^^Lu-DOTATATE SPECT 
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Purpose Patient-specific dosimetry of lutetium-177 
(^^^Lu)-DOTATATE treatment in neuroendocrine tumours 
is important, because uptake differs across patients. 
Single plioton emission computer tomography 
(SPECT)-based dosimetry requires a conversion factor 
between tlie obtained counts and tlie activity, wliicli 
depends on tlie collimator type, the utilized energy 
windows and the applied scatter correction techniques. 
In this study, energy window subtraction-based scatter 
correction methods are compared experimentally and 
quantitatively. 

Materials and methods ^^^Lu SPECT images of a 
phantom with known activity concentration ratio between 
the uniform background and filled hollow spheres were 
acquired for three different collimators: low-energy high 
resolution (LEHR), low-energy general purpose (LEGP) 
and medium-energy general purpose (MEGP). Counts 
were collected in several energy windows, and scatter 
correction was performed by applying different methods 
such as effective scatter source estimation (ESSE), 
triple-energy and dual-energy window, double-photopeak 
window and downscatter correction. The intensity ratio 
between the spheres and the background was measured 
and corrected for the partial volume effect and used to 
compare the performance of the methods. 



Results Low-energy collimators combined with 208 keV 
energy windows give rise to artefacts. For the 113keV 
energy window, large differences were observed in the 
ratios for the spheres. For MEGP collimators with the ESSE 
correction technique, the measured ratio was close to the 
real ratio, and the differences between spheres were small. 

Conclusion For quantitative ^^^Lu imaging MEGP 
collimators are advised. Both energy peaks can be utilized 
when the ESSE correction technique is applied. The 
difference between the calculated and the real ratio is less 
than 1 0% for both energy windows. NucI Med Commun 
35:522-533 © 2014 Wolters Kluwer Health | Lippincott 
Williams & Wilkins. 
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Introduction 

One of the treatment options for neuroendocrine 
tumours is peptide receptor radionuclide therapy with 
lutetium-177 (^^^Lu) -DOTATATE, which contains the 
octreotide derivative (Tyr'^)-octreotate in its molecular 
structure. It will bind to the somatostatin receptors, 
which overexpress in these types of tumours [1]. The 
combination of |3-radiation and y-radiation makes ^^^Lu 
an interesting choice for this treatment. |3-Radiation is 
necessary for therapy, whereas y-radiation makes imaging 
of the therapy (uptake) possible. Ideally, a dose as high as 
possible is given to the tumour, whereas the surrounding 
organs receive a dose as low as possible. In ^^^Lu- 
octreotide treatments, the uptake of the tracer differs 
across patients, which influences the dose to the 
neighbouring organs as well. Therefore, in order to obtain 
optimal treatment, patient-specific dosimetry is necessary. 
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Knowledge about dose response in peptide receptor 
radionuclide therapy is still limited, mainly because there 
continue to be considerable uncertainties in dose 
calculation [2]. The method used to perform lutetium- 
octreotide dosimetry varies considerably [3-5], which 
makes it difficult to compare treatment outcome and 
complication rates depending on the administered dose. 
In order to perform accurate dosimetry a conversion 
factor is necessary between the obtained counts in a 
single photon emission computer tomography (SPECT) 
image and the real activity. 

^^^Lu has two main photon energy peaks that could be 
used for imaging and therefore for dosimetry: one at 
208keVand one at 113keV Low-energy collimators are 
normally suitable only for energies below 160keV When 
isotopes with a higher photon energy are used, septal 
penetration will occur, which can give rise to star 
artefacts. This makes low-energy collimators unsuitable 
for imaging the 208 keV photons of ^^^Lu. Because the 
208 keV peak has a higher incidence (10.4%) compared 
with that of 113keV (6.8%) and in addition a slightly 
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Fig. 1 




The 10-1 acylic phantom (dimensions 24.1 x 30.5 x 24.1 cm) with six 
spheres of inner diameters 1 0, 1 3, 1 7, 22, 28 and 37 mm. The phantom 
is described in NEMA NU2-2007, section 7 'Image quality accuracy 
of attenuation and scatter corrections'. Image courtesy of PTW 
(Physikalisch-Technische Werkstatten). 



lower attenuation in the tissue, it is more suitable for 
imaging, and therefore normally medium-energy general 
purpose (MEGP) collimators are used for imaging with 
^^^Lu. Many centres use only the higher energy peak 
because of downscatter from the 208 keV peak into the 
113keV window. This can originate from scatter in the 
patient and backscatter into the crystal, as 180° scattering 
of 208 keV photons gives rise to 1 15 keV photons. Further, 
at Rigshospitalet, at present only the high-energy peak is 
utilized for imaging, combined with an MEGP collimator 
and the standard scatter correction [effective scatter 
source estimation (ESSE)]. 

Acquiring data from the 113keV window in addition to 
the 208 keV window will improve the count statistics and 
thus might also increase contrast resolution, which makes 
more accurate dosimetry or shorter scanning time 
possible. Appropriate scatter correction techniques have 
to be applied in order to be able to use both windows. 

Several methods exist for scatter correction - for 
example, energy window subtraction methods [6-11], 
convolution subtraction methods [12] and limitation of 
the detection of scattered photons by, for example, 
changing the energy window [7]. In this study different 
(down) scatter correction techniques, relying on energy 
window subtraction, were compared with each other for 
three different collimators in an experimental ^'^'^Lu 
SPECT study to find the optimal situation for dosimetry 
in ^^^Lu imaging. In a study by Beauregard et al. [13] one 
energy window subtraction method with medium-energy 
collimators was investigated for the 208 keV window. It is 
hypothesized that quantitatively accurate ^'^'^Lu SPECT 
can be achieved using both energy windows when the 
appropriate combination of collimator and scatter correc- 
tion technique is used. 



Materials and methods 

A 10-1 acrylic phantom (NEMA2007/IEC2008) with six 
fillable spheres (inner diameter: 37, 28, 22, 17, 13 and 
10mm; wall thickness: ^ 1 mm) was filled with a ^'^'^Lu 
solution (Fig. 1). The hollow spheres were filled with a 
higher activity concentration than the uniform back- 
ground activity concentration, resulting in hot spots. Two 
samples were removed from the background and from the 
largest of the spheres to perform accurate counting in a 
calibrated gamma counter (Packard Cobra II Auto 
Gamma; Perkin-Elmer, Waltham, Massachusetts, USA). 
The true intensity ratio between the spheres and the 
background was calculated on the basis of these 
measurements. The absolute activity concentration in 
the spheres at the beginning of the first scan was 
1.31MBq/ml, and the sphere-to-background concentra- 
tion ratio was 12.88. These activity concentrations are 
comparable to the activity concentration in normal 
clinical practice; a homogeneous distribution of 7.4 GBq 
^^^Lu in a 70 kg patient gives a concentration of 
^O.lMBq/ml, which was used as the background 
concentration in the phantom. Tumours can have an 
activity concentration of the order 1 MBq/ml or higher, 
and kidneys typically have an uptake on the order of 
0.5 MBq/ml 1 day after administration. 

Experiments were performed using a dual-head 16-slice 
computer tomography (CT) Philips Precedence SPECT/ 
16MDCT scanner (Philips Healthcare, Best, the 
Netherlands) with a 9.5-mm-thick scintillation crystal. 
SPECT scans were consecutively acquired within 3 h with 
three different collimators (in this order): low-energy 
high resolution (LEHR), MEGP and low-energy general 
purpose (LEGP). Table 1 shows the specifications of the 
collimators. Scans were acquired over 360° at 128 angles 
(30s/angle, step and shoot mode) in a 128 x 128 matrix 
size with a pixel size of 3.195 mm and an orbit radius 
of 25 cm. The time interval between the starting times of 
the three scans was ^ 1 h. Fifteen different energy 
windows were defined (Fig. 2 and Table 2). Two main 
energy windows were defined, one at 113keVand one at 
208 key both with a 20% width. In order to obtain the 
appropriate energy windows some of the acquired energy 
windows in Table 2 need to be added and/or subtracted 
before applying the energy window-based correction 
methods described in the next section; for example, the 
upper half of the photo peak at 208 keV can be obtained 
by subtracting window 15 from window 1. Reconstruction 
for all scans and all collimators was performed with the 
Philips Astonish reconstruction technique. This is a 
technique based on ordered subset expectation max- 
imization with a depth-dependent resolution recovery 
method [14]. Reconstruction was performed with three 
iterations and eight subsets as recommended by the 
vendor. In order to investigate the possible influence 
of the number of iterations and subsets, some of the 
reconstructions for the LEHR collimator were also 
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Table 1 Collimator specifications according to product sheet 



Collimator names 


Hole size 


Hole length 


Septa thickness 


Relative sensitivity^ 


System resolution FWHM at lOcm"^ 


Septal penetration (%)'^ 


LEHR 


1.40 


32.8 


0.152 


0.55 


7.4 


1.3 


LEGP 


1.40 


24.7 


0.180 


1.0^ 


8.8 


1.9 


MEGP 


2.95 


48.0 


1.143 


0.8^ 


11.3 


3.3 



All collimators have hexagonal hole shape. All sizes in mm. 

FWHM, full width at half maximum; LEHR, low-energy high resolution; LEGP, low-energy general purpose; MEGP, medium-energy general purpose. 
^Compared with LEGP for ^^'^Tc, 20% energy windows. 
'^For a 9.5-mm-thick crystal. 

^Septal penetration of respectively 140keV photons for the low-energy collimators and of 300 keV photons for the medium-energy collimators. 
'^265cpm/^LCi = 1 19cps/MBq. 
^Measured ratio for ^^""Tc. 




Table 2 The acquired energy windows and their description 



Window numbers 


Lower energy limit (keV) 


Upper energy limit (keV) 


Description 


Abbreviation 


1 


187.6 


229.2 


Main window 1 




2 


101.7 


124.3 


Main window 2 




3 


177.6 


187.6 


TEW 1 


TEWwide 


4 


229.2 


239.2 


TEW 1 


TEW^ide 


5 


183.6 


187.6 


TEW 1 


TEWsmaii 


6 


229.2 


233.2 


TEW 1 


TEWs^all 


7 


91.7 


101.7 


TEW 2 


TEW^ide 


8 


124.3 


134.3 


TEW 2 


TEW^ide 


9 


97.7 


101.7 


TEW 2 


TEWsmall 


10 


124.3 


128.3 


TEW 2 


TEWe^all 


11 


146.0 


187.6 


Scatter 1 


DE\N_k factor 


12 


79.1 


101.7 


Scatter 2 


DE\N_k factor 


13 


208.4 


250.0 


Downscatter 1 


DS 


14 


124.3 


146.9 


Downscatter 2 


DS 


15 


187.6 


208.4 


Dual peak 1 


DPW 



DEW, dual-energy window; DPW, double-photopeak window; DS, downscatter; TEW, triple-energy window. 



performed with more than the default three iterations 
and eight subsets - namely, with three iterations and 
16 subsets and with four iterations and 16 subsets. 
One of the reconstructions was performed with eight 
iterations and 16 subsets. 



Neither prefiltering nor postfiltering was applied. At- 
tenuation correction was applied using a low-dose 
140 kVp GT scan [15], which was transformed into 
an attenuation coefficient map for 100 keV photons, 
which subsequently was scaled to an attenuation map for 
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the 113 and 208 keV photons. As scatter correction 
techniques are appHed, narrow beam attenuation coeffi- 
cients are appHed. 

Scatter correction methods 

In this section the different scatter correction methods 
are described using the following symbols: 

(1) 6" for absolute number of counts; 

(2) w for the window width in keV; 

(3) c for the count concentration in counts/ml; 

(4) V for the volume in ml; 

(5) k for the multiplication factor. 

The subscripts prim and total are used to refer to the 
primary (unscattered) photons and the total number of 
counts, respectively. 

The following scatter correction methods were used in 
this study. 

No scatter correction 

For comparison we started with a basic SPECT recon- 
struction in which only attenuation correction was 
applied, based on the acquired GTscan [15]. No scatter 
correction was applied. 

Effective scatter source estimation 

The standard Philips scatter correction, which is based on 
the ESSE method, was also applied [16]. This method 
uses the concept of defining an effective source distribu- 
tion of scatter, calculating the expected scatter projection 
data from this effective source of scatter photons and 
adding this calculated scatter projection data to the 
projection of the estimated image to produce the total 
estimated projection. The estimate of the scatter is 
calculated from the image estimate and density map by 
convolving the image estimate with a three-dimensional 
scatter convolution kernel. This convolution kernel is 
precalculated and depends on the photon energy [16]. 
No extra energy window, besides the main energy 
window, is needed for this correction method. 

Triple-energy window 

With the triple-energy window (TEW) method two small 
energy windows are defined around the main energy peak 
window. The two abutting windows are used to estimate 
the amount of scatter in the main window. Scatter 
correction was applied following [8,17]: 

^prim — ^total 7^ (Slower ^ Slower ~l~'^higher ^ ^higher) 5 

(1) 

K\ — , 

where the subscript lower refers to the lower energy 
window and the subscript higher refers to the higher 
energy window. The subscript i is used to indicate an 
arbitrary window. 



In the literature, different energy window widths are 
used [8,18-20]. Therefore, in this study two different 
versions of the TEW method were investigated. For this 
purpose two different sets of energy windows were 
defined, one with a width of 10 keV (windows 3 and 4, 
and windows 7 and 8) and one with a width of 4keV 
(windows 5 and 6, and windows 9 and 10). 

Duai-energy window 

In the dual-energy window (DEW) method a broad 
window is defined below the main energy window with 
the same width. The assumption is that the spatial 
distribution of the scattered photons in the lower energy 
window is a good estimate of the distribution of the 
scattered photons in the main energy window. Subtrac- 
tion is then applied in the following way [6]: 

^prim — ^total^-^s X^S5 (2) 

where the subscript s is used to indicate the scatter 
window. 

In the literature there is some discussion on the value 
of k that should be used. In the original paper of Jaszczak 
<?/ a/. [6] a factor of 0.5 was used for ^^"^Tc, but in later 
papers, for example, in those by Koral et al. [21] and 
Luo et al. [ZZ]/it was stated that the multiplication factor 
is dependent on, among other things, the ratio of 
background activity (sphere vs. cylinder) and the volume 
of interest (VOI) size. Monte Carlo simulations sug- 
gested a multiplication factor of 0.43 [23]. The multi- 
plication factors found in these studies were all for ^^"^Tc. 
In this study we used three different multiplication 
factors to determine their influence on the corrected 
image: 0.5, 0.8 and 0.9. The value of 0.5 was used because 
it is the value most used in the literature, and the values 
of 0.8 and 0.9 were used on the basis of the ratio of the 
background activity and the VOI sizes as suggested by 
other studies such as the ones by Koral and collea- 
gues [21,22,24]. 

Downscatter correction 

To correct for downscatter a window above the main 
energy window was defined with the same width as the 
main energy window (windows 13 and 14). The corrected 
image was calculated as follows [25]: 

^prim — ^total ~^ds ^ ^ds 7 (^) 

where the subscript ds is used to indicate the down- 
scatter window. 

To determine ^^s a VOI outside the phantom was drawn. 
In the projection data outside the phantom no primary 
photons are present. In contrast, downscattered photons 
can be detected if they penetrate through the septa of 
the collimator, while it is assumed that they are 
collimated perpendicular to the crystal surface. Alterna- 
tively, they can be backscattered into the crystal 
and detected outside the projection of the phantom. 
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Assuming that k^^ is space invariant (not depending on 
position), was calculated using the formula: 

E ^total 



In the case of imaging the 113 keV photons of ^'^'^Lu, the 
208 keV photons can downscatter into the 113keV 
window. The 208 keV photons can scatter both inside 
the phantom and in the material behind the scintillation 
crystal. Backscattering (180°) of a 208 keV photon results 
in a photon of 115 key which is well inside the imaging 
window for the 113keV photons. 

Downscatter correction for the few percent 250 and 
321 keV photons is not necessary in the case of imaging 
the 208 keV photons with medium-energy collimators 

Combined scatter and downscatter correction 

With the TEW method, correction for scatter and 
downscatter is applied, but with the DEW method 
only a part of the downscatter is taken into account. 
To compensate for downscatter as well, a combination 
of the DEW method and the downscatter correc- 
tion method was applied. The used windows were 
146-187.6 keV (window 11) with a i'-factor (kz) of 0.5 
for the scatter correction and 208.4-250 keV (window 2) 
for the downscatter correction with the same ^-factor 
(>fi) as used for downscatter alone for the 208.4 keV 
photon energy peak and 79.1-101.7 (window 12) and 
124.3-146.9 (window 14) for the 1 13 keV photon energy 
peak. Because the scatter window is contaminated with 
downscatter as well, a correction factor ki, is applied to 
correct for this. This correction factor ^3 is calculated in 
the same way as the i'-factor for downscatter [see Eq. 
(4)], except for the fact that instead of the main energy 
window 6'pj-im the scatter window Cs was used [26]. The 
final formula for obtaining the corrected image then 
becomes: 



prim 



C'totai - {h -h xh)Cds -h^C^- 



(5) 



Dual-pliotopeak window 

Instead of measuring extra windows next to the main 
energy window to estimate the scatter in the main 
window, the information in the main window itself can be 
used. The dual-photopeak window (DPW) method is 
based on the fact that the lower part of the main energy 
window contains more scattered photons than the upper 
part (ratio a). In the study by Ljungberg et a/. [27] the 
main energy window is split into two abutting nonover- 
lapping energy windows that are symmetrically located 
around the photopeak. The observed relative difference 
between the two windows is then used to calculate the 
scatter fraction, or the scatter-to-total ratio (STR). In this 
study we have used a method that is comparable to the 
one used in the study by Pretorius era/. [28] and is called 
the channel ratio method. Assuming that the contribution 



of the nonscattered photons to the upper and lower 
window is the same {Cpnm/2), it can be derived from 



prim 



^total ~^total STR, 



with 



orr-iT-, 0C~I~1 Slower ^higher 
STR = 7 X ^ 



(6) 



(7) 



a— 1 Slower H" ^higher 

The value of a (the ratio between the scatter in the lower 
and the upper window) equals 3 in the case of a linear 
descending distribution of the scatter in the primary 
window. 

The method by Pretorius et al. [28] is based on 
experimental determinations of the window weights, 
whereas Ljungberg ^/ [27] have calculated a fit (based 
on Monte Carlo simulations) from which a can be 
determined. We determined a based on comparison with 
the formula used by Ljungberg et al. [27] and it was 
calculated to be 3.6, which is close to the value of a linear 
descending distribution. 

Because the dual-peak method is highly noise sensi- 
tive [9], a spatial 3x3 filter with a Gaussian kernel 
(coefficients 1/4, 1/8 and 1/16) is applied before 
calculating the scatter fraction. 

Because only 15 energy windows could be defined in the 
Philips scanner and because we do not expect the DPW 
to work properly in the 113 keV window because of the 
downscatter of the 208 keV photons that violates 
the assumption of the scatter distribution, we have 
chosen to define the DPW window only for the 208 keV 
window. 

Analysis 

vols were drawn around the six spheres and in the 
background of the phantom (Fig. 3). 

In order to compensate for the partial volume effect and 
subjectivity, large volumes were drawn as proposed by 
Tossici-Bolt al. [29]. The counts in the spheres 
(^sphere) wcrc Calculated with: 

^sphere = ^VOI -^Bg(^VOI- ^sphere ) , (8) 

where the subscript BG refers to the background. Note 
that the volume of the sphere is the outer volume 
including the cold walls. For all the different correction 
methods, the intensity ratios between the six spheres and 
the background were calculated and compared with the 
real value. The number of counts in the background 
outside the phantom was also obtained. If the correction 
methods work correctly, there should not be any counts 
there. 

To be able to calculate the absorbed dose in a patient a 
conversion factor between the number of counts in the 
organ of interest and the activity is needed. This 
conversion factor is calculated in a phantom in which 
the activity and the uptake (100%) are known. To 
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Fig. 3 




Two transaxial slices (medium-energy general purpose collimators and 208 keV energy window), reconstructed and both attenuation and scatter 
corrected with Astonish, (a) Low-dose computer tomography slice through the centre of the (visible) walls of the spheres and the volume of interest 
(vols) in different colours, (b) Single photon emission computer tomography (SPECT) of the same slice as in (a) with the VOIs shown in different 
colours, (c) SPECT slice in the phantom far away from the spheres illustrating the VOI for the background. The banana-like VOI outside the phantom 
is used for determination of the weight for certain energy window correction methods. Colour scales are linear with highest intensity at the top. 



calculate the conversion factor the same VOIs are used as 
for calculating the ratios. The real activity was calculated 
on the basis of the measurements in a gamma counter, 
which in turn is calibrated against the dose calibrator 
used for the measurement of injected activity in 
phantoms and patients. Activity values are corrected for 
the decay time and the efficiency of the gamma counter 
itself. The conversion factor is obtained by dividing the 
real activity by the measured concentration in the sphere, 
corrected for the background activity. 

Statistics 

In order to compare the different methods with the 
reference value (the measured intensity ratio) analysis of 
variance tests were performed. If any significant {P < 0.05) 
difference was found between the groups a post-hoc 
Dunnett test (two-sided) was performed to compare with 
the reference value. Bonferroni tests were performed to 
compare the different variations of the methods with each 
other (e.g. TEW with small and broad windows). When the 
data were not normally distributed (tested with 
Kolomogorov-Smirnov) the Kruskal-Wallis test was applied 
instead of the analysis of variance. When the variances in 
the data differed significantly from each other, the Dunnett 
T3 test was performed instead of the Dunnett test. 

Results 

Comparison between the different scatter correction 
methods was performed for each collimator separately. 
Because the inaccuracy of the intensity ratio becomes 
larger when the sphere becomes smaller (because of 
partial volume effects and the limited amount of counts), 
the smallest two spheres are not taken into account for 
the comparisons. Figure 4 shows three SPECT reconstruc- 
tions for both the 113 keV window with LEGP collimators 
and the 208 keV window with MEGP collimators to 
illustrate the influence of scatter and downscatter 
correction. Note that downscatter correction is not 
relevant for the 208 keV window. Figure 5 shows three 



SPECT reconstructions for both the 113keVand 208 keV 
energy windows in order to illustrate the differences 
between the three collimators. 

The results for the ratios are shown in Figs 6-9. Because 
of the high values for the ratios for DEW0.9 these results 
are omitted from the figures, but are shown in Tables 3 
and 4. 

No significant differences were found between the two 
TEW methods (one with a small and one with a somewhat 
broader window) for any of the three collimators {P=\ .000) . 

For all the correction methods, the ratio between the 
counts outside the phantom and the background inside the 
phantom (a measure for the contribution of downscattered 
photons to the final images) was smaller than 10%, except 
for the TEWsmaii method in combination with the LEHR 
(12%) for the 113 keV window. For the MEGP collimator 
this ratio was below 1% for all correction methods. 

There were no significant differences between the ratios 
for the images reconstructed with four and eight instead 
of three iterations and with 16 instead of eight subsets. 

Significant differences were found between the DEW 
methods with different ^-factors. Depending on the 
collimator used, this may result in overestimation and 
underestimation of the ratio. 

Low-energy high resolution 

The reconstruction of the higher energy window (around 
208 keV) with the LEHR collimator gives artefacts, which 
results in ratios more than a factor 2 lower than the real 
value. Scatter correction methods were not applied to this 
window because of this problem. 

Significant differences between the calculated ratio and the 
real ratio were found for the 208 keV window and for the 
113keV window without scatter correction. High relative 
SDs in the background volume were observed, with values 
above 15% even rising up to 49% for the DEW method with 
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Fig. 4 




Six single photon emission computer tomography reconstructions [low-energy general purpose collimators for 1 13keV (a, b, and c) and medium- 
energy general purpose for 208 keV (d, e, and f)] illustrating the effect of scatter and downscatter correction, (a, d) Slice through the centre of the 
spheres with attenuation correction only, (b, e) Same as (a) and (d) but with scatter correction with the dual-energy window method (weight factor of 
0.5). (c, f) Same as (a) and (d) but with combined scatter and downscatter correction with the triple-energy window method with 'wide' energy 
windows. Colour scales are linear with highest intensity at the top. 



Fig. 5 




Reconstructions of the three collimators [low-energy high resolution (a, d), low-energy general purpose (b, e) and medium-energy general purpose 
(c, f)] for 113 keV (a, b, c) and 208 keV (d, e, f) energy windows. All slices are reconstructed with attenuation correction and the effective scatter 
source estimation method. Colour scales are linear with highest intensity at the top. Notice the low image quality for the 208 keV energy window with 
low-energy collimators due to septal penetration. 



a >f-factor of 0.9. For the LEHR collimator and the 113keV 
imaging window the DEW0.5 method and the combined 
downscatter and DEW method give similar values, because 



the downscatter correction in the DEW window is close to 
the amount of downscatter in the 113keV imaging 
window. Figure 6 shows an overview of the results. 
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Fig. 6 



20.00 
18.00 



Sphere 1 □ Sphere 2 ■ Sphere 3 □ Sphere 4 




Results of the low-energy high resolution collimator for 113 keV. DEW, dual-energy window; DS, downscatter; ESSE, effective scatter source 
estimation; TEW, triple-energy window. 



Fig. 7 



Sphere 1 □ Sphere 2 ■ Sphere 3 □ Sphere 4 



20.00 
18.00 
16.00 




Results of the low-energy general purpose collimator for 1 13keV. DEW, dual-energy window; DS, downscatter; ESSE, effective scatter source 
estimation; TEW, triple-energy window. 



Low-energy general purpose 

As for the LEHR collimator the reconstruction of the 
208 keV window does not work well because of septal 
penetration of the primary imaging photons. Reconstruc- 
tion with the scatter correction techniques was not 
applied because of this. For the LEGP collimator and 
the 113 keV imaging window the DEW0.5 method and the 
combined downscatter and DEW method give exactly 



the same values, because the downscatter correction in 
the DEW window equals the amount of downscatter 
in the 113 keV imaging window. 

The ratios obtained with 208 keV differ significantly from the 
real ratio. For the 113 keV window, only the method without 
scatter correction gives significantly different results from the 
real value. Figure 7 shows an overview of the results. 
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Fig. 8 
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Results of the medium-energy general purpose collimator for 1 1 3 keV. DEW, dual-energy window; DS, downscatter; ESSE, effective scatter source 
estimation; TEW, triple-energy window. 
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Results of the medium-energy general purpose collimator for 208 keV. DEW, dual-energy window; DS, downscatter; ESSE, effective scatter source 
estimation; TEW, triple-energy window. 



Medium-energy general purpose 

When the MEGP collimator is used the calculated ratios 
for the higher energy window are already close to the real 
ratio without applying scatter correction. When the 
standard scatter correction (ESSE) is used the calculated 
ratio for both the lower and the higher energy window 
is close to the real ratio. The other scatter correction 
methods tend to overestimate the ratio, especially in 



the lower energy window. The downscatter method does 
not work in the higher energy window, as there is no 
downscatter present in this situation, and it does 
not work in the case of the downscatter correction for 
113 keV as in the study by de Nijs al. [25], as the 
high-energy photons are collimated, and the weight 
cannot be determined outside the patient. Relative 
SDs of the background are high, especially for the 
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Table 3 Difference (%) between tlie real intensity ratio and the ratio of the largest sphere estimated after the different scatter correction 
methods for the three collimators 





LEHR 


LEHR 


MEGP 


MEGP 


MEGP 


LEGP 


LEGP 


Energy window (keV) 


113 


208 


113 


208 


113 + 208 


113 


208 


None 


-42.5±9.1 


-66.0±13.2 


-35.0±7.2 


-10.9±9.7 


-271 ±7.1 


-44.1175 


-72.1±13.7 


ESSE 


-18.6±11.7 


-64.7±12.5 


-8.9±9.0 


-6.6±9.4 


-7.8±9.1 


-21.3±7.4 


-71.2±13.2 


TEWsmall 


-12.9±25.8 




94.2±40.9 


10.5±18.5 




-4.3±23.7 




TEWwide 


-5.4±29.7 




107.0±46.1 


12.0±16.6 




3.7±31.7 




DEW, k=0.5 


-18.4±17.0 




53.8±26.3 


15.0±17.7 




-14.9±15.2 




DEW, k=0.8 


21.7±35.9 




180.5±61.8 


40.9±26.6 




46.9±48.5 




DEW, k=0.9 


47.3±47.4 




285.4±97.1 


52.7±29.2 




94.2±75.8 




DS1 


-23.0±16.1 




66.3±69.0 


-102.0±73.6 




-19.9±12.6 




DS1 +DEW0.5 


-16.6±19.7 




47.9±22.7 


16.3±16.7 




-14.9±15.2 




DPW 








-8.0±10.5 









The statistical uncertainty due to the noise of all calculated ratios is less than 2 percentage points. A minus sign indicates an underestimation. Standard deviation of the 
differences for the four largest spheres is given as ±SD. 

DEW, dual-energy window; DPW, dual-photopeak window; DS, downscatter; ESSE, effective scatter source estimation; LEHR, low-energy high resolution; LEGP, 
low-energy general purpose; MEGP, medium-energy general purpose; TEW, triple-energy window. 



Table 4 Conversion factor, in kBq/cps, of the largest sphere estimated after the different scatter correction methods for the three 
collimators 





LEHR 


LEHR 


MEGP 


MEGP 


MEGP 


LEGP 


LEGP 


Energy window (keV) 


113 


208 


113 


208 


113-^208 


113 


208 


None 


7.2±1.7 


6.8±12.6 


5.1 ±0.7 


7.6±0.9 


3.1 ±0.4 


4.5±0.8 


4.5±8.4 


ESSE 


11.3±2.2 


8.1 ±9.7 


8.2±1.0 


8.9±1.1 


4.2±0.5 


7.1 ±0.8 


5.3±6.2 


TEWsmall 


12.0±9.5 




12.1 ±4.6 


8.9±2.4 




7.2±3.9 




TEWwide 


11.2±10.9 




11.2±4.7 


8.8±2.0 




6.7±6.6 




DEW, k=0.5 


9.3±3.5 




9.2±2.5 


8.5±2.0 




5.8±1.7 




DEW, k=0.8 


12.0±10.4 




11.7±5.0 


9.4±3.1 




7.4±8.5 




DEW, k=O.Q 


13.5±15.0 




13.1 ±75 


9.8±3.2 




8.6±18.3 




Dsr 


9.3±3.3 










5.7±1.3 




DS1 +DEW0.5 


9.7±4.7 




8.9±2.1 


7.8±1.7 




5.8±1.7 




DPW 








7.4±1.1 









The statistical uncertainty due to the noise of all calculated conversion factors is less than 2 percentage points. Standard deviation of the conversion factor for the four 

largest spheres is given as ± SD. 

DEW, dual-energy window; DPW, dual-photopeak window; DS, downscatter; ESSE, effective scatter source estimation; LEHR, low-energy high resolution; LEGP, 
low-energy general purpose; MEGP, medium-energy general purpose; TEW, triple-energy window. 
^The downscatter correction method does not work for collimated high-energy photons (MEGP). 



113keV window. Figures 8 and 9 show an overview of 
the results. 

Table 3 shows an overview of the difference between the 
real intensity ratio and the estimated intensity ratio 
after the different scatter correction methods are applied 
for the three collimators used in this study. 

Table 4 shows an overview of the conversion factors 
needed to calculate the observed activity concentration 
in a patient from the number of counts in the organ of 
interest after the different scatter correction methods are 
applied for the three collimators used in this study. 

Equation (8) is used to calculate the statistical un- 
certainty due to the noise, expressed as the empirical SD 
of the calculated intensity ratio, based on VOI statistics. 
For the VOI statistics the data from the final recon- 
structed images are used, which include uncertainties 
due to, for example, counting statistics and reconstruc- 
tion. For spheres 1-4 (which are the four largest spheres) 
this statistical uncertainty due to noise typically spanned 
from less than 0.1-0.4, respectively, except for the ratios 
determined with the TEW method for the MEGP 



collimator with 113keV windows, in which these values 
were 0.2 and 0.7, respectively The uncertainty for the 
DEW method with 113keV windows was up to 1.8 for 
sphere 4 because of failure of this method. However, the 
relative uncertainty was still less than 8%. The statistical 
uncertainty due to noise in spheres 5 and 6 can be more 
than 20% depending on the correction method. Image 
deterioration caused by noise was not investigated. 

Discussion 

In this study different scatter correction methods based 
on energy window subtraction are compared with each 
other in an experimental ^^^Lu SPECT study for three 
different collimators. For ^''''Lu imaging the MEGP 
collimator is advised. Both energy peaks could be used 
when the standard scatter correction technique (ESSE) 
was applied. The difference between the calculated ratio 
and the real ratio was less than 10% for both the low- 
energy and the high-energy window and when both 
windows were combined. Using one appropriate conver- 
sion factor (common for all spheres) both the 208 and the 
113keV imaging window can be utilized on their own or 
combined after reconstruction. 
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As expected, the reconstruction of the higher energy 
window (208 keV) is incorrect for the two low-energy colH- 
mators because of septal penetration. This leads to 
'leaking' of information into neighbouring pixels, which 
renders accurate calculation of ratios and doses impos- 
sible. Applying scatter correction techniques is useless in 
this case. 

Large differences were observed in the ratios calculated 
for the different spheres, especially for the low-energy 
window (113keV). The high relative SDs in the back- 
ground influence the accuracy of the obtained concentra- 
tions. Small differences in the background concentration 
can highly influence the obtained concentration ratio, 
especially in the smaller spheres, which might explain the 
differences between the spheres. Differences in the ratio 
of the spheres also mean different conversion factors for 
different sizes of spheres. In patients, the size of the 
tumour is normally not known, which implies that only 
one conversion factor, suitable for all tumour/organ sizes, 
can be used. Therefore, only methods that yield the same 
value for each sphere can be used for accurate dosimetry. 

For some of the methods, several variations were also 
tested. For the TEW method, two different energy 
windows widths were tested, one smaller with a width of 
4keVand one somewhat broader, with a width of lOkeV 
Both methods performed comparably. Because the TEW 
method is noise sensitive [7], which is especially a 
problem in dynamic SPECT studies [30,31] in which few 
counts are acquired, the TEW method with the broader 
windows might be more stable for dynamic studies. 

For the DEW method different i'-factors were tested. 
Statistically significantly different results were obtained 
with these methods. Which method performs best 
depends on which collimator is used. From the litera- 
ture [21,22] it has become clear that the i'-factor used 
depends on different factors. Koral et al. [21] found that 
the i'-factor decreased as the background concentration 
increased. The size of the VOI also influences the 
i'-factor, but its influence is dependent on the background 
concentration. When the background concentration is low, 
the k value drops rapidly with increasing VOI size, but k is 
almost constant when the background concentration is 
large. With the VOI size and background-to-sphere 
activity that were used in this study, a ^-factor of around 
0.5 should be appropriate. Unfortunately, Koral et al. [21] 
investigated the i'-factor only for ^^"^Tc and an LEGP 
collimator. Luo et al. [ZZ] investigated the ^-factor for 
99mrp^ and ^^^I, but a different phantom was used, which 
influences the value of the i'-factor. Our study of ^^^Lu 
suggests that for the 113keV photons a ^-factor of 0.5 for 
all collimators might be feasible. This i'-factor of 0.5 also 
holds for the 208 keV photons with the MEG? collimator. 
The disadvantage of optimizing this i'-factor is that this is 
only possible with experimental or simulation studies and 



that the values obtained might differ from the values that 
are needed in patients because the ratio between the VOI 
and the background in patients is not known beforehand. 

Standard (down) scatter correction methods such as TEW 
and DEW do not work very well for lutetium for various 
reasons. The problem with the DEW method in the 
113keV window might be that characteristic radiation 
from the collimator lead (energies between 72 and 
88 keV) deteriorates the DEW In the case of imaging of 
the 208 keV photons the TEW method might be changed 
to a DEW method with one small window, because the 
upper window will mostly contain noise, which influences 
the accuracy of the method. Downscatter correction in 
the 208 keV window is not really necessary because 
almost no downscatter will be present in this window as a 
result of the very low abundance (few per cent compared 
with the 208 keV photons) of the higher energy photons. 
When correction with this method is applied, the 
calculated ratio differs considerably from the real ratio. 
This might be due to large uncertainties arising from the 
noise in determining the >f-factor when no downscatter is 
present. When the MEGP collimator is used, large 
differences occur between the different scatter correction 
methods, but most of them are not significantly different 
from the real value. This might be due to the large 
variances obtained. 

The DPW for the 208 keV window and the standard 
scatter correction technique (ESSE) for both the 113 and 
the 208 keV window perform best with MEGP collima- 
tors. When ESSE is used, the difference between the real 
ratio and the calculated ratio is small, and the differences 
between the spheres are small. When the DPW method is 
used, the differences between the different sized spheres 
are larger, which is reflected by a larger SD. 

Only 15 energy windows could be defined in the Philips 
scanner. Because we did not expect the DPW to work 
properly in the 113keV window as a result of the down- 
scatter of the 208 keV photons that violates the assumption 
of the scatter distribution, we chose to define the DPW 
window for the 208 keV window only. The DPW method 
is sensitive to deviations of the energy calibration of the 
scanner [32], as it is based on dividing the main energy 
window into two parts with the limit placed exactly at the 
imaged photon energy, and it is noise sensitive [9]. In 
practice, this can be an important source of error. 
Furthermore, the DPW method is possibly unstable [33]. 

A ratio for the smaller spheres (mainly spheres 5 and 6) 
cannot be obtained, which means that dosimetry in small 
tumours/organs will fail as well. There are two main 
reasons for this: the influence of the partial volume effect 
will be larger in smaller objects because of the resolution 
limit of the scanner, and the number of counts will 
be lower, which makes the measurement less precise. 
The phantom corresponds to a patient of about 70 kg. 
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As the resolution of the reconstructed images depends 
on the colHmator source distance, the resolution will be 
worse for larger patients and consequently the minimal 
detectable tumour size will be larger. 

For ^^^Lu imaging the MEGP collimator is advised. Both 
113 and 208 keV energy windows can be used when the 
standard scatter correction technique (ESSE) is applied. 
The difference between the calculated ratio and the real 
ratio is less than 10% for both energy windows and when 
both windows are combined after reconstruction. When 
this method is used one conversion factor can be used for 
all the spheres, which makes accurate dosimetry feasible. 
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